Phase diagram of 4 He on graphene 



Jodok Happacher 

Department of Physics, Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, 
University of Munich, Theresienstrasse 37, 80333 Munich, Germany 

Philippe Corboz 
Theoretische Physik, ETH Zurich, CH-8093 Zurich, Switzerland 



o 

O 

Q 



O 

o 



> 

ON 

(N 



X 



Massimo Boninsegni 
Department of Physics, University of Alberta, Edmonton, Alberta, Canada, T6G 2G7 

Lode Pollet 

Department of Physics, Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, 
University of Munich, Theresienstrasse 37, 80333 Munich, Germany [j 
(Dated: January 1, 2013) 

The low temperature phase diagram of 4 He adsorbed on a single graphene sheet is studied by 
computer simulation of a system comprising nearly thousand helium atoms. In the first layer, two 
commensurate solid phases are observed, with fillings 1/3 and 7/16 respectively, separated by a 
domain wall phase, as well as an incommensurate crystal at higher coverage. No evidence of a 
thermodynamically stable superfliuid phase is found for the first adlayer. Second layer promotion 
occurs at a coverage of 0.111(4) A -2 . In the second layer two phases are observed, namely a 
superfluid and an incommensurate solid, with no commensurate solid intervening between these two 
phases. The computed phase diagram closely resembles that predicted for helium on graphite. 



I. INTRODUCTION 

The study of films of highly quant al fluids, such as he- 
lium, is motivated by the search for novel phases of mat- 
ter in confinement and/or reduced dimensions. Indeed, 
an experimentally controllable way of making quasi-2d 
interacting 4 He systems is to adsorb a thin film of 4 He 
on a substrate. On weakly attractive substrates, such as 
those of some alkali metals, 4 He forms superfluid films 
(down to a monolayer 1 thin) whose thickness smoothly 
increases with chemical potential, with no evidence of 
layering^ On the other hand, on stronger substrates 
adsorption occurs through the formation of successive, 
well-defined layers (up to seven on graphite), with essen- 
tially no quantum-mechanical atomic exchanges taking 
place between the first few adjacent adlayers; the phase 
diagram on these substrates is richer, displaying a variety 
of phases, including crystalline ones, either commensu- 
rate or incommensurate with the underlying substrate. 

Crowell and Reppy 4 raised some time ago the possibil- 
ity of a "supersolid" phase, - characterized by simultane- 
ous density and superfluid long-range order, in the vicin- 
ity of a possible crystalline phase of the second helium ad- 
layer, registered with the underlying graphite substrate. 
This contention has recently been reiterated^ The most 
reliable, first-principle numerical studies of helium films 
on graphite have yielded no evidence of such a phase, as 
no registered crystal is observed in the second adlayer. 7 

Graphene (a single sheet of graphite) has also been 
theoretically considered as a possible substrate for helium 
adsorption, and the phase diagram of the adsorbate in the 
low temperature (i.e., T — » 0) limit has been computed 
by means of different numerical techniques^— A single 



sheet of carbon atoms is somewhat less attractive than 
a graphite substrate; quantitatively, the atomic binding 
energy for the first 4 He adlayer is reduced by approxi- 
mately 10% (about 13.4 K) compared to graphite^ One 
might imagine that the energy offset could lead to dif- 
ferent physical behavior, but first principle calculations 
suggest that a difference of that order of magnitude in the 
adsorption potential is likely to have little or no effect on 
the phase diagram of helium on graphite £ Thus, one may 
expect no quantitative differences in the phase diagram 
of 4 He adsorbed on graphene versus that on graphite. 
Indeed, that is the conclusion at which most numerical 
studies carried out so far for this system have arrived, 
with an outstanding puzzle concerning a possible super- 
fluid response in the first adsorbed layer, near and at 
commensurate filling, reported in Refs. ^ and [III. 

We report in this article results of a theoretical study 
of the low temperature phase diagram of 4 He adsorbed 
on graphene, based on computer simulations. We con- 
sider here a system comprising at least twice as many 
4 He atoms than in previous studies by others, the goal 
being that of attempting a reliable extrapolation of the 
physics of the system in the thermodynamic limit. The 
main findings of our study are largely in line with most 
previous works on graphene, but with no evidence of any 
"supersolid" phase, either in the second nor in the first 
adsorbed layer. 

The remainder of this paper is organized as follows: 
first, we briefly discuss the model of the system that is 
used, as well as the computational methodology; then, 
we proceed with the illustration of the results, focusing 
on the first and second adlayers. We summarize our main 
results in the Conclusions, where we also address issues 
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that may be the subject of future work. 



II. METHODOLOGY 

In order to study numerically the physical proper- 
ties of 4 He on a graphene sheet, we performed equi- 
librium, large-scale computer simulations of a model of 
the system of interest, using the continuous-space Worm 
Algorithm i 12 i 13 Graphene is modeled as an ideal, two- 
dimensional honeycomb lattice, with a carbon-carbon 
bond length a=1.42 A. Carbon atoms are treated as fixed 
particles in our simulation, an assumption justified by 
their relatively large mass, compared to that of the 4 He 
atoms. The system is enclosed in a three-dimensional 
cell, shaped as a parallelepiped; the graphene sheet is 
aligned parallel to the xy plane (at z=0). Periodic 
boundary conditions are used in all directions, but the 
box is sufficiently elongated in the z direction to make the 
boundary condition immaterial. The ensuing, quantum- 
mechanical many-body Hamiltonian is the following: 



i<j iL 



where m is the mass of the 4 He atoms, v is a pair- 
wise potential describing the interaction between two he- 
lium atoms, whereas v describes the interaction of each 
helium-carbon pair, with the position of a carbon 
atom. Both u and v are assumed to depend only on 
the relative distance between two particles. The interac- 
tion between helium atoms v is taken to be the accepted 
Aziz pair potential^ whereas for the carbon-helium po- 
tential u a Lennard- Jones potential is used, with pa- 
rameters e and a chosen following Ref. [l5|. The simu- 
lated systems comprise close to one thousand 4 He atoms. 
The largest honeycomb lattice simulated here has size 
49.19 Ax 51.12 A, and consists of 960 carbon atoms. All 
the simulations are performed in the grand canonical en- 
semble, at finite temperatures ranging from 0.5 K to 1.0 
K; the results at the two different temperature differ very 
little, suggesting that they are essentially ground state es- 
timates, at least as far as energetics and structures are 
concerned. Finally, all the results presented here are in- 
dependent of the initial atomic configuration utilized in 
the simulation. 



III. RESULTS FOR THE FIRST LAYER 

We start the discussion of our results with the first 
layer. Fig JT] shows the computed average equilibrium 4 He 
coverage (two-dimensional density) as a function of chem- 
ical potential fi. 

As the chemical potential is increased above /i = — 130 
K, the average 4 He density first jumps from zero to the 
value 0.0636 A -2 corresponding to the crystalline, com- 
mensurate Cl/3 phase, wherein one of three equivalent 
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FIG. 1. (Color online). Average coverage versus chemical 
potential for the first 4 He adlayer on graphene. Data shown 
are for different simulations, corresponding to different ini- 
tial atomic arrangements. The first layer becomes populated 
around \i — —130 K, forming a Cl/3 crystal phase. The hor- 
izontal lines denote the density corresponding to maximum 
first layer occupation, the C7/16 phase and the Cl/3 phase, 
top to bottom. 



adsorption sites of the lattice is occupied by a single 4 He 
atom, as seen in Fig. [2j It is worth mentioning that in 
some simulations a metastable liquid phase with only a 
slightly higher value of the energy per atom was observed; 
conceivably, such a phase might be seen in experiments. 16 
No evidence of a low-coverage superfluid liquid (which ex- 
ists on a smooth, flat substrate 1 ^) is observed. 

The Cl/3 crystal forms at a distance z = 2.864(7) A, 
with a layer width Az = 0.388(9) A(these numbers hold 
for the entire first layer), only slightly greater than that 
for graphite. The thermodynamic stability of the Cl/3 
phase is signaled by an extended plateau of the density 
as a function of /i, hinting at a large gap in the spec- 
trum with zero superfluid response and zero compress- 
ibility. All of this is very similar to what is observed 
on graphite. 7 A direct calculation of the superfluid frac- 
tion, supplemented with the observation that particle ex- 
changes are strongly suppressed in the crystalline phase, 
clearly shows that its value is zero in the thermodynamic 
limit, for the Cl/3 commensurate phase, which is con- 
sistent with what is also reported in Ref. [Ill, as well as 
with results of lattice model calculations 

Ref. 9 reports a very small (less than 1%), but fi- 
nite value of the superfluid fraction at T=0 in corre- 
spondence of the Cl/3 crystalline phase, increasing to 
approximately 14% upon doping such phase with vacan- 
cies. Because that of Ref. [9| is a ground state calculation, 
a direct comparison with our results is not possible, as 
the transition temperature of a hypothetical superfluid 
crystalline phase could plausibly be much lower than 0.5 
K, which is the lowest temperature considered here. Two 
remarks are in order, however: a) it is far from clear to 
us whether the computational methodology adopted in 
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Ref. 0, namely Diffusion Monte Carlo, indeed allows one 
to resolve reliably a value of the superfluid fraction less 
than 1%; the authors do not discuss at all its robustness 
against major sources of systematic errors (finite system 
size and finite population bias^), either of which could 
easily account for the reported superfluid response, and 
b) the thermodynamic stability of commensurate phases 
doped with vacancies is also not discussed in Refs. |9| and 
[ill . It is worth mentioning that grand canonical simu- 
lations (carried out in this work) clearly have an edge 
over those performed in the canonical ensemble, when it 
comes to ruling out phase coexistence. Our simulations 
are consistent with the general observation that any crys- 
talline phase doped with either vacancies or interstitials 
is metastable in this kind of system. 

Further increasing the chemical potential leads to the 
appearance of domain walls, as is seen in Fig. [2j which 
first occur along one principal axis of the Cl/3 phase. 
This phase is akin to striped phases, but with increas- 
ing chemical potential a proliferation of domain walls is 
observed, along more and more principal axes. In the 
thermodynamic limit, it has no superfluid response. This 
process stops when a commensurate crystalline phase la- 
belled as C7/16 with a coverage of 0.0835 A -2 is formed, 
and which can be seen as every unit cell of the Cl/3 phase 
being surrounded by domain walls. The C7/16 structure 
has a rhombic unit cell, with 7 4 He atoms distributed 
over 16 absorption sites. For this commensurate struc- 
ture, too, no finite superfluid signal is computed. Fig. [2] 
shows the 4 He density profile for the C7/16 phase ob- 
served in our large scale simulations. Surprisingly, the 
helium atoms are in our case distributed more uniformly 
over the unit cell than in previous works iii Moreover, 
the structure is slightly rotated (a comparison with the 
commensurate phases seen by others is shown in Fig.[3)). 
We performed the simulation with a smaller system size 
(the same as the one used in Ref. [Tlh , and we also used as 
initial configuration one corresponding to the the density 
profile obtained in Ref. [U, but our Monte Carlo simula- 
tion still stabilized, after a sufficiently long time, to the 
same structure of higher symmetry inside the unit cell, 
shown in Fig. [3j Thus, we conclude that the difference 
is not attributable to system size, nor is it a result of a 
specific choice of initial atomic configuration. 

The most likely reason for the difference between our 
and previous results is the modeling of the helium-carbon 
interaction. While we use here a two-body potential (u 
in Eq. [1]) with fixed carbon positions, thereby retain- 
ing full rotational symmetry, previous studies made use 
of an external one-body potential obtained by summing 
anisotropic Lennard- Jones potentials proposed to fit he- 
lium scattering from graphite. The slightly different ways 
of accounting for corrugation may well be at the root of 
the lower symmetry structure found in previous studies, 
and this is supported by comparing the actual potentials 
felt at locations inside the unit cell, which differ between 
the two methods. 

On increasing /i even further, the C7/16 phase is re- 




FIG. 2. (Color online). Top Row: Left: The commensurate 
Cl/3 phase. Top Row: Right: Domain wall structure at a 
coverage 0.072 A -2 . Bottom Row: Left: C7/16 commensurate 
phase as found in this work (see also Fig. [3]). Bottom Row: 
Right: 4 He density profile for the incommensurate crystalline 
phase, at a coverage 0.111 A -2 , which is the maximum first 
layer coverage 




FIG. 3. (Color online). Left: 4 He density profile for the 
C7/16 phase as seen in Ref. [Til . The unit cells are shown 
with dashed lines. The small (white) dots represent the un- 
derlying graphene lattice. The large (red) dots are maxima 
of the density. Right: 4 He density profile for the C7/16 phase 
obtained in this study. 



placed by an incommensurate, compressible solid phase, 
sketched in Fig. [2] (panel on the left in the bottom row). 
The density increases linearly with chemical potential un- 
til a density of 0.107(3) A -2 is reached, where second 
layer promotion begins. This value is slightly higher than 
the one reported in Ref. [U, but lower than the one from 
Gordillo et al, namely 0.115 A -2 —. It is also lower than 
the corresponding one for a graphite substrate, either 
from theoretical calculations, 7 namely 0.114 A -2 , or in- 



- 


, , ♦ t * 


- 


i ^ * 




1 

1 










-26 


-24 -22 -20 -18 -16 -14 




|l(K) 



FIG. 4. (Color online). Average coverage versus chemical 
potential for the second 4 He adlayer on graphene (offset by the 
maximum first layer density, m = 0.111 A -2 ). Data shown 
are for different simulations, corresponding to different initial 
atomic arrangements. The second layer becomes populated 
around ji = — 25 K, forming a homogeneoous superfluid. The 
horizontal line denotes the maximum density in the first layer 
before layer promotion starts; this first layer density was also 
seen when the second layer is strongly occupied. For higher 
//, a first-order phase transition to an incommensurate solid 
occures. No evidence of any "supersolid" phase is found. 

ferred from specific heat measurements^ i.e., 0.127 A -2 , 
something that can be attributed to the weaker substrate 
attraction. 



IV. RESULTS FOR THE SECOND LAYER 

We now turn to the second layer. The first stable phase 
that we observe is a uniform liquid, which is superfluid 
at the temperature considered here (Fig. [4] and Fig. [5]). 
The equilibrium coverage is roug hly 0.150 A" 2 , consis- 
tent with the two-dimensional density of the superfluid 
being close to the equilibrium density of two-dimensional 
4 He, i.e. 0.043 A -2 . This superfluid phase has a rela- 
tively wide domain of existence, extending up to a layer 
density of 0.076 ± 0.003 A -2 . Within the quoted statis- 
tical uncertainty, this value is in agreement with that of 
Ref. [l(| and also consistent with that found on graphite^ 
namely 0.076 A -2 . 

On raising the chemical potential, a first-order phase 
transition to an incommensurate crystal is observed, 
which persists up to atomic third layer promotion. The 
density profile of the incommensurate solid phase is 
shown in Fig. [5] Its superfluid response is zero in the 
thermodynamic limit. This is qualitatively the same be- 
havior observed for graphite (see Ref. 7). Nowhere in 
the second layer are commensurate phases found, and in 
particular we find no evidence of any "supersolid" phase. 
The arguments for the absence of the suggested C7/12 
and C4/7 phases formulated in Ref. 7 for graphite apply 
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FIG. 5. (Color online). Left: second layer density profile 
at a two-dimensional density 0.0731 A -2 . The dark spots 
correspond to the underlying helium first layer atoms. Right: 
second layer density profile for the incommensurate phase, at 
a two-dimensional density 0.0855 A -2 , which is immediately 
below third layer promotion. 



equally well for graphene. Gordillo and Boronat sug- 
gested that there could be a coexistence region between 
liquid and incommensurate phase. 10 Since the transition 
between a liquid and solid where the substrate does not 
play a role has to be first order, such a region must be 
sought around fi = — 18K in Fig. |4j We did not an- 
alyze this first order transition in more detail, since it 
was very narrow already in graphite 1 As a final remark, 
we note that the 4 He density between the layers is es- 
sentially zero, i.e., layers constitute distinct, effectively 
two-dimensional systems. 



V. CONCLUSION 

We studied the phase diagram of the first and second 
layer of 4 He adsorbed on graphene. In the first layer, 
we find a very stable Cl/3 phase, followed by a domain 
wall phase, a C7/16 phase and an incommensurate solid. 
This is the same as for graphite^ and in line with previ- 
ous studies jSill The only difference we observed consisted 
of the density distribution of the helium atoms inside the 
C7/16 unit cell, where we observed a small rotation (and 
more symmetric) distribution. We attribute this differ- 
ence to the modeling of the carbon-helium interaction, 
which was done slightly differently than in previous stud- 
ies: We used a fully rotationally symmetric two-body 
potential between carbon and helium atoms, and kept 
a (classical) hexagonal carbon lattice in the simulation, 
while previous studies used a one-body potential that 
integrated out the C atoms. We neglected defects, dis- 
order, and other lattice imperfections in the graphene. 
For the second layer, our results are also in line with the 
ones found on graphite 7 and previous studies. We find 
first the liquid-gas coexistence, then a superfluid and an 
incommensurate solid. No commensurate structures are 
found, nor any supersolids. Layer promotion occurs at 
a slightly lower chemical potential than for graphite, be- 
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cause of the slightly weaker helium-substrate potential. 

In future work, we will study more realistic descrip- 
tions of graphene, which is not entirely flat but has intrin- 
sic microscopic roughening. Suspended graphene sheets 
under a transmission electron microscope were seen to 
have a surface normal varying by several degrees and the 
out of plane deformations reach 1 nm^ This roughened 
structure is thought to suppress and stabilize thermal vi- 
brations and to influence the behavior as an adsorbent. 



Our setup can also be modified to study helium inside 
carbon nano-tubes: instead of using periodic boundary 
conditions in both x and y directions on the graphene 
substrate, one could use periodic boundary conditions 
only in the x direction. 
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